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SIMULATION  AND  PERFORMANCE  OF  BRUSHLESS  DC  MOTOR  ACTUATORS 

SUMMARY 

The  simulation  model  for  a  Brushless  D.C.  Motor  and  the  associated 
commutation  power  conditioner  transistor  model  are  presented.  The  necessary 
conditions  for  maximum  power  output  while  operating  at  steady-state  speed  and 
sinusoidally  distributed  air-gap  flux  are  developed. 

Comparisons  of  simulated  model  with  the  measured  performance  of  a  typical 
motor  are  done  both  on  time  response  waveforms  and  on  average  performance 
characteristics.  These  preliminary  results  indicate  good  agreement.  Plans 
for  model  improvement  and  testing  of  a  motor-driven  positioning  device  for 
model  evaluation  are  outlined. 
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_  LIST  OF  SYMBOLS 

Bm  Viscous  friction  coefficient  of  the  motor  (oz-in/rad  per  sec) 

B]_  Viscous  friction  coefficient  of  the  load  (oz-in/rad  per  sec) 

EBA  Back  EMF  of  phase  A  (volts) 

EBB  Back  EMF  of  phase  B  (volts) 

EBC  Back  EMF  of  phase  C  (volts) 

IA  Phase  A  current  (amperes) 

IB  Phase  B  current  (amperes) 

Iq  Phase  C  current  (amperes) 

Im  Power  supply  current  (amperes) 

J^  Moment  of  inertia  of  motor  shaft  (oz-in-sec2) 

Jl  Moment  of  inertia  of  load  shaft  (oz-in-sec2) 

kb  Back  EMF  constant  (volts/rad  per  sec) 

kt  Torque  constant  (oz-in/ampere) 

L  Inductance  of  the  stator  winding  (Henrys) 
IAB,  IqA,  Igc   Loop  currents  (amperes) 
kadi  Air-gap  flux  adjustment  factor 

kwl  Magnitude  of  speed  versus  torque  curve  slope  (rad  per  sec/oz-in) 

Ra  Resistance  of  the  stator  windings  (ohms) 

Rs  Power  supply  interval  resistance  (ohms) 
<I>A»  ^B'  4>C   Per  Pnase  air-gap  flux  (webers) 

Rq  'On'  resistance  of  the  power  transistor  (Ohms) 

Tm  Motor  restraining  torque  (oz-in) 

T;l  Load  restraining  torque  (oz-in) 
TA,  TB,  Tq  Per  phase  developed  torque  (oz-in) 

wm  Speed  of  the  motor  shaft  (rad/sec) 


GR   Gear  ratio 

9    Angular  displacement  (radians) 

VA,  Vg,  Vq  Phase  terminal  voltage  (volts) 

V5   Power  supply  terminal  voltage  (volts) 

RPSA,  RPSB,  RPSC   Position  sensors  logic  signals 

PWM  Pulse  width  modulation 


INTRODUCTION 

Recent  improvements  in  rare-earth  magnetic  materials  for  use  in  brushless 
dc  motors  have  allowed  reconsideration  of  electro-mechanical  actuator  systems 
for  applications  requiring  very  high  ratios  of  torque-to-inertia.  The 
investigation  discussed  herein  has  been  concerned  with  characterizing  mathe- 
matically the  dynamical  features  of  a  missile  fin  actuation  system,  from  the 
input  to  the  brushless  dc  motor  to  the  output  shaft  of  the  mechanical  actuator. 
The  physical  model  is  based  upon  an  existing  prototype  actuator  currently 
under  evaluation  at  the  Naval  Weapons  Center,  China  Lake,  California. 

In  general,  brushless  dc  motors  produce  torque  through  the  interaction  of 

a  magnetic  field  generated  by  a  permanent  magnet  rotor  and  a  dc  generated 

magnetic  field  in  the  stator.  The  rotating  permanent  magnet  eliminates  the 

rotating  armature  and  the  mechanical  wear  normally  associated  with  brushes. 

These  motors  fall  in  the  class  of  Permanent  Magnet  Motors  and  enjoy  certain 

advantages  over  wound-field  types  such  as: 

"...Linear  torque-speed  characteristics,  high  stall  (accelerating) 
torque,  no  need  for  electric  power  to  generate  the  magnetic  flux  and 
a  smaller  frame  and  lighter  motor  for  a  given  output  power"  [1]. 

Additionally,  the  brushless  dc  motor  is  characterized  by: 

"...controllability  over  a  wide  range  of  speeds,  capable  of  rapid 
acceleration  and  deceleration,  convenient  control  of  shaft  speed  and 
position,  no  mechanical  wear  problem  due  to  commutation  and  better  heat 
dissipation  arrangement"  [1]. 

The  fundamental  requirement  of  an  electro-mechanical  actuator  control 

system  is  to  provide  torque  to  an  output  shaft,  sense  the  position  of  the 

shaft  and  adjust  the  torque  to  balance  the  load  when  the  desired  position  is 

reached.  This  must  be  accomplished  with  a  minimum  of  frictional  resistance 

and  delays  associated  with  the  inertia  of  the  mechanical  components.  Effects 

of  viscous,  static  and  coulomb  friction,  together  with  the  torque  required  to 

accelerate  the  mechanical  components  of  the  system,  lead  to  a  reduction  in 


torque  available  at  the  output  shaft  and  an  associated  reduction  in  system 
performance. 

One  approach  to  the  analysis  of  the  electro-mechanical  actuator  system 
has  been  to  divide  the  system  into  two  sequential  problem  areas.  The  first 
deals  with  the  dynamic  analysis  of  the  brushless  dc  motor  and  development  of 
the  transfer  function  necessary  to  duplicate  actual  steady-state  and  transient 
performance.  The  second  area  deals  with  modeling  the  mechanical  system 
elements,  taking  as  input  the  dc  motor  shaft  angular  acceleration  predicted  by 
the  motor  analysis.  The  mechanical  system  must  be  modeled  considering  the 
effects  of  friction  and  inertia  and  translating  the  rotational  motion  of  the 
brushless  dc  motor  shaft  to  the  output  shaft  of  the  actuator  for  application 
to  missile  maneuvering  control.  The  results  obtained  from  a  study  that  has 
placed  primary  emphasis  upon  the  latter  problem  area  -  the  modeling  of  the 
mechanical  drive-train  leading  to  the  fin  shaft  was  presented  in  Ref.  3  and  4. 

This  report  deals  with  the  first  problem  area  and  presents  an  overview  of 
the  work  done  by  MacMillan  [4].  In  Ref.  4,  MacMillan  developed  a  model  for 
the  output  circuitry  of  a  transistorized  power  conditioner  that  provided  the 
required  motor  commutation.  The  model  developed  in  Ref.  4  used  an  ideal  (zero 
impedance)  power  supply.  This  report  extends  the  model  to  include  internal 
power  supply  resistance  and  also  identifies  areas  for  further  improvement  in 
the  system  modeling. 

The  next  section  of  this  report  presents  a  brief  description  of  the 
system  followed  by  a  development  of  the  model  and  an  analysis  of  the  simulated 
motor  performance.  Results  of  the  simulated  system  are  then  compared  to  the 
measured  response  from  a  typical  commercially  available  motor  and  recommen- 
dations for  further  improvements  in  the  model  are  outlined. 


.SYSTEM  DESCRIPTION 
GENERAL 

The  system  is  viewed  as  a  position  control  device  to  maintain  an  output 
angle  under  an  applied  hinge  moment  due  to  aerodynamic  forces  on  a  fin  or 
aileron.  The  motor  is  a  permanent  magnet  dc  motor  with  feedback  in  the  form 
of  back  emf  proportional  to  the  angular  velocity  of  the  motor.  The  block 
diagram  of  the  dc  motor  and  power  conditioner  is  shown  in  Fig.  1  [4].  The 
mechanical  actuator  and  drive  train,  as  currently  envisioned,  introduce 
various  inertial  and  damping  loads  together  with  an  aerodynamic  force  and  its 
associated  fin  hinge  moment  that  must  be  overcome  to  produce  output  motion. 
An  operational  block  diagram  of  the  load  torque  is  shown  in  Fig.  2  where  the 
hinge  moment  and  motor  shaft  angular  acceleration  are  viewed  as  inputs  to  the 
drive  train  [2].  Figure  3  is  a  schematic  of  the  drive  train  which,  as 
presently  constituted,  includes  the  motor  shaft  (leadscrew),  ball  screw 
assembly,  and  the  crank  which  is  keyed  to  the  output  (fin)  shaft.  Inertial 
loads  are  considered  individually  within  three  major  subdivisions  of  the 
actuator;  the  output  shaft  to  crank,  crank  to  ball  screw,  and  ball  screw  to 
leadscrew. 

A  detailed  analysis  and  development  of  the  mechanical  model  with  the 
associated  assumptions  are  contained  in  Reference  2.  As  stated  in  the 
introduction,  this  report  reviews  the  development  of  the  motor  model  reported 
in  Reference  4.  For  convenience  in  the  development  of  the  model,  MacMillan 
considered  the  load  seen  by  the  motor  to  be  caused  by  the  inertia  and  kinetic 
(viscous)  friction  of  the  motor  shaft  as  described  by  Thomas  in  Reference  5. 

The  schematic  diagram  of  the  power  supply,  power  conditioner  and  motor- 
load  is  shown  in  Figure  4  where  it  is  noted  that  the  power  conditioner  logic 
and  driver  circuits  are  modeled  as  ideal  on-off  devices  with  zero  time  delay. 
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Figure  2     Load  Torque  Block  D 


lagram 


Gp-LL  screw 


Figure  3     Sketch  of  Actuator  Arrangement 
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Figure  4   Schematic  Diagram  of  Power  Supply,  Power  Conditioner 
and  Motor  with  Inertia  and  Viscous  Friction  Load. 
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Additionally,  MacMillan  considered  the  power  supply  to  be  ideal  (zero  interval 
resistance)  to  further  reduce  the  complexity  and  simplify  the  development  of 
the  power  conditioner  model.  However,  in  this  report,  the  model  used  for  the 
power  supply  includes  the  internal  resistance,  Rs.  The  effect  of  including  Rs 
in  the  requirements  and  performance  of  the  system  model  is  presented  in  the 
section  on  Analysis  of  Simulated  Motor  Performance.  For  additional  informa- 
tion on  the  development  of  the  D.C.   motor  model  refer  to  Reference  4. 

DEVELOPMENT   OF   SYSTEM  MODEL 

The  basic  simulation  by  Thomas  used  a  single  power  supply  and  super- 
imposed the  phase  currents  to  produce  the  motor  torque.  Motor  drive  was  then 
realized  by  multiplying  averaged  armature  current  by  a  torque  factor. 

In  the  model  developed  by  MacMillan,  the  supply  is  considered  to  be  a 
split  supply  of  equal  voltages.  Armature  current  is  assigned  a  positive  sign 
if  it  flows  in  the  positive  direction  (i.e.  into  the  motor).  The  use  of  a 
split  power  supply  allowed  MacMillan  to  apply  circuit  reduction  techniques. 
Because  of  the  symmetry  that  resulted  from  the  split  power  supply  approach, 
the  complex  3~phase  bridge-type  circuit  simplified  into  a  two-window  network. 
The  development  of  the  reduced  circuit  begins  with  the  definitions  of  loop 
currents  as  shown  in  Figure  5  where  the  power  supply  voltage  of  2V  is 
represented  as  V+  and  V-.  The  simulated  model  equivalent  circuit  is  shown  in 
Fig.  6  with  node  N  defined  as  the  mid-point  of  the  power  supply  and  no*de  0  to 
be  the  mid-point  of  the  motor  windings  and  neither  of  these  nodes  are 
considered  to  be  at  ground  potential  (the  ground  is  identified  as  the 
Reference  node).  Note  that  the  power  supply,  VIN  =  VIF  -  VIB  -  2V  and  that 
all  the  network  variables  are  defined  using  the  CSMP  model  variable 
definitions   [6].      Given  the  assumption  that  the  network   is    balanced,    that    is, 
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all  transistors,  diodes  and  field  windings  are  approximated  as  being  alike  in 
characteristics,  it  is  possible  to  apply  the  circuit  reduction  technique  of 
Thevenin  and  obtain  the  simple  2-window  network  shown  in  Figure  7  [4],  In 
Figure  7,  the  Thevenin  equivalent  voltages  and  currents  are  defined  in  terms 
of  the  CSMP  model  variables. 

The  back  emf  voltages  of  each  phase  (VEMFA,  VEMFB  and  VEMFC)  are 
developed  and  presented  in  Reference  H  in  which  it  is  shown  that  these 
voltages  are  summed  2  at  a  time  in  proper  sequence  to  compute  the  loop 
currents  and  resulting  phase  torque.  The  total  developed  torque  is  then  the 
sum  of  the  2-active  winding  phase  torques  over  the  proper  60  degree  of 
mechanical  angle.  Figure  8  shows  the  back  emf  voltages  generated  across  2 
windings  taken  2  at  a  time.  The  60  degree  of  mechanical  angle  that  produces 
torque  is  shown  by  the  logic  unit  output  levels  of  the  position  transducer 
RPSA,  RPSB  and  RPSC  in  Figure  8. 

Reference  4  presents  a  detailed  development  of  the  power  conditioner 
model  that  includes  the  assumptions  for  the  switching  transistor  dynamics  used 
in  the  power  conditioner  model  as  well  as  the  development  of  the  harmonic 
air-gap  flux  used  in  the  motor  model. 

ANALYSIS  OF  SIMULATED  MOTOR  PERFORMANCE 
The  proper  evaluation  of  a  system  model  requires,  quite  naturally,  a 
comparison  of  the  actual  system  response  for  the  same  type  of  input  function 
used  for  the  model.   For  the  brushless  D.C.  motor,  it  is  assumed  that 
measurements  can  be  obtained  for  the  following  voltages  and  currents: 

-  Power  supply  terminal  voltage;  VSg 

-  Voltage  across  two  of  the  three  windings;  Va^ ,  V^q ,  and  Vca 
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Figure  6     Two-Window  Equivalent  Circuit  of  the  Balanced  Bridge  Network 
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-  Voltage  from  one  winding  terminal  to  ground;  Va>  Vb  and  Vc 

-  Power  supply  current;  Im 

-  Phase  current;  Ima,  Imb,  Imc 

Before  taking  data  from  the  simulated  model,  it  is  necessary  to  set  all 
constants  and  parameters  to  the  nominal  values  given  by  the  supplier  of  the 
motor.  Care  must  be  taken  to  properly  interpret  the  values  given  by  the 
manufacturer. 

In  particular,  the  value  of  the  Back  EMF  coefficient,  kb,  needs  special 
attention.  If  the  motor  data  sheet  indicates  that  this  coefficient  is 
obtained  by  measurement  of  peak  to  peak  voltage  across  two  windings,  then  in 
the  model  that  uses  distributed  sinusoidal  air-gap  flux,  an  adjustment  factor, 
kadi  ^s  required  to  insure  that  the  Back  EMF  values  are  in  agreement.  It 
follows  that  the  torque  coefficient,  k^,  must  also  be  adjusted  by  this  same 
scale  factor  [7].  The  value  of  ka(^j  can  be  computed  from  the  measured  no-load 
data  for  the  motor  current,  speed  and  applied  voltage  by  using  the  average 
characteristics  balance  equation, 


lm  ■  ^vsg 


kbwnl)/(Ra  +  Rq) 


where  kb  =  kbm  kadj 

kbm  =  manufacturer   supplied  Back  EMF  constant. 
Thus   the  adjustment   factor  becomes 

kadj    ■    l>sg  -   im(Ra   +   Rq)]/(kbmwnl) 

For    example,    given  a  no-load  speed   of   3060   RPM  with  a  terminal   voltage  of 

30  volts  and  a  no-load  current  of  0.30   amperes,    kacjj    =    0.825    (given   a   total 

series   resistance   of    1.47   ohms).      The   value   of    ka<jM   must  of  course  be  less 

than   unity,    otherwise   the  motor   current   would   be   negative    in   value   which 

implies   that  generator   rather   than  motor  action   is   taking  place. 
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Before  an  attempt  is  made  to  close  the  loop  on  the  motor-load  to  form  a 
positioning  device,  it  is  important  to  operate  the  system  as  a  velocity  device 
and  obtain  measurement  of  the  motor  average  performance  characteristics. 
Typical  motor  average  performance  curves  are  speed  versus  torque,  motor 
current  versus  torque  and  output  power  versus  torque  from  no-load  torque  to 
near  stall  torque  conditions.  Figure  9  shows  these  three  curves  for  the  model 
with  torque  given  in  oz-in  units.  Note  in  particular  the  straight  line 
characteristics  for  motor  speed  and  current  and  that  the  power  output  is  of 
quadratic  form.  The  condition  of  constant  input  voltage  is  used  when 
gathering  data  for  these  curves,  and  it  follows  that  if  speed  is  a  straight- 
line  versus  torque,  then  current  also  has  a  straight-line  when  plotted  against 
torque.     This  can  be  shown  as  follows: 

wm  "  wnl  ~  kwlT 
where  wm  =  motor  speed 

wnl  =  no-load  speed 

kwl  =  magnitude  of  slope  of  the  speed  vs  torque  curve 
T  =  load  torque 
The  current,  using  average  values,  can  be  written  as 

im  =  <vsg  "  kbwm)/(Ra  +  V 
where  vsg  =  power  supply  terminal  voltage 

kb  =  Back  EMF  coefficient 

Ra  =  winding  resistance 

Rq  =  transistor  "ON"  resistance 
then  by  substitution, 

im  =  [(vsg  -  kb  wnl)/(Ra  +  Rq)]  +  [(kbkwl)/(Ra  +  Rq)]T 
where  the  first  term  on  the  right  hand  side  of  the  equation  represents  the 
no-load  current. 
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Figure  9  Average  Performance  Characteristics  of  the  Simulated  Model 
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The  power  output  versys  torque  curve  in  Figure  9  reaches  a  peak  value 
between  no-load  torque  and  stall  torque.  The  analytical  development  for  this 
quadratic  curve  is  presented  in  Appendix  A  where  it  is  shown  that  the 
condition  for  peak  power  output,    P0,    is  that  the  torque  must  be  of  value 

T   =    (2kbwnl   -  vsg)/(2kbkwl) 
to  produce  maximum  power  output 

po  =  ktvsg2/(4Rskb) 

Since   the   ratio  k^/kb   is  a  constant,    increasing  peak  power  output  depends 

upon  decreasing  the  power  supply  interval  resistance,  Rs,  or  increasing  the 
power  supply  voltage,  vsg  as  one  would  logically  expect.  It  must  be 
understood  that  the  assumption  used  in  the  above  analysis  was  that  the  air- 
gap  flux  was  a  constant  average  value.  For  the  actual  motor  the  air-gap  flux 
is  distributed  sinusoidally  and  current  in  the  windings  produces  a  field  that 
may  cause  distortion  in  the  air-gap  flux.  At  large  current  values  that  occur 
at  and  above  peak  power  output,  the  field  distortion  may  result  in  an  increase 
or  a  reduction  of  output  power  for  a  given  load  torque. 

Another  motor  coefficient  to  consider  carefully  is  the  no-load  viscous 
friction  constant,  Bm .  The  value  of  Bm  is  often  not  included  in  the  motor 
data  sheet  and  the  probable  reason  is  that  its  value  will  depend  somewhat  on 
the  manner  in  which  the  motor  is  attached  to  the  system  load.  Bm  can  be 
calculated  from  given  no-load  data  as  follows: 

Bm  -   <kt  kadj    im)/wnl 
where  im  is  the  no-load  current. 

RESULTS  AND   CONCLUSIONS 
The    preliminary    results    indicate   that    the   balanced    bridge   circuit 
approach  used   in  the  development  of  the  Brushless  DC  motor  model  produces   good 
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agreement  with  measured  motor  characteristics  as  indicated  below.  Further 
evaluations  of  the  model  as  well  as  improvements  and  additions  to  the  model 
will  be  conducted  in  the  near  future  using  data  currently  being  gathered  from 
a  prototype  of  a  fin  positioning  actuator. 

A  comparison  of  model  and  motor  produced  Back  EMF  voltage  is  shown  in 
Figure  10.  The  upper  curve  is  the  simulated  motor  waveform  and  the  lower 
curve  is  the  measured  Back  EMF.  The  model  adjustment  factor,  kacj j  was  set  to 
0.63  value  by  trial  and  error  until  peak  to  peak  voltages  were  in  agreement. 
Another  verification  of  the  model  is  obtained  by  comparison  of  Figure  11  and 
Figure  8.  Both  figures  show  waveforms  of  Back  EMF  across  2  windings  taken  2 
at  a  time  and  also  show  the  timing  waveforms  for  the  Position  Sensor  Devices 
(RPSA,  RPSB  and  RPSC)  for  counter  clockwise  rotation.  The  waveforms  agree  in 
both  phasing  and  in  form. 

The  steady-state  performance  curves  of  a  typical  motor  are  given  in 
Figure  12  where  the  motor  load  torque  is  given  in  lb  in  units.  These  curves 
were  produced  with  a  constant  terminal  voltage  of  30  volts  and  agree  in  form 
with  the  curves  of  the  model  as  given  in  Figure  9.  Peak  power  output  occurs 
at  load  torque  of  approximately  95  oz-in  for  both  motor  and  model. 

Additional  validation  of  the  model  is  shown  in  Figure  13  where  the 
current  in  Phase  C  for  the  model  (upper  curve)  is  in  close  agreement  both  in 
form  and  in  phase  with  the  same  current  for  the  motor  (lower  curve). 
Additional  model  results  are  shown  in  Appendix  B  where  typical  input  voltage 
step  response  waveforms  are  presented. 

FURTHER  RESEARCH 
Improvements  in  the  model  are  required  with  regard  to  the  triggering  on 
of  the  diodes  used  to  protect  the  transistors  from  excessive  reverse  currents 
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Figure  10     Back  EMF  Voltage:     Top  Curve-Model    Output, 
Bottom  Curve-Typical   Motor  Output 
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Figure  12     Typical   Measured  Motor  Performance 
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Figure  13     Phace  C  Current:      Top  Curve-Simulated  Mddel, 
Lower  Curve-Typical   Motor 
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during  the  switch  off  time  for  the  transistors.  MacMillan  in  Ref.  4  was  able 
to  provide  this  action  for  the  condition  that  the  power  supply  interval 
resistance,  Rs  is  zero  value.  When  Rs  is  not  zero,  additional  programming  is 
required  to  perform  the  required  trigger  action. 

The  Pulse  Width  Modulation  (PWM)  used  by  Askinas  in  Ref.  8  must  be  added 
to  the  present  model,  however  Askinas  used  a  separate  additional  transistor  to 
perform  PWM.  The  Power  Conditioner  to  be  used  in  the  Test  Stand  for 
measurement  of  actuator  performance  accomplishes  PWM  by  using  the  lower  set  of 
the  commutating  power  transistors.  Therefore  some  modification  of  Askinas' 
CSMP  coding  may  be  required   to  accurately  perform  the  PWM  speed  control. 

Extensions  of  the  model  will  include  a  position  sensing  device  and  a 
tachometer  as  well  as  the  nonlinear  loading  as  developed  by  Franklin  in  Ref  9. 
Test  data  results  will  be  compared  to  model  in  both  transient  (step)  response 
and  in  frequency  (Bode)  response  and  the  results  will  be  reported  in  a  Summary 
Report  for  FY1986. 
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APPENDIX  A 
Power  Output  of  a  Brushless  DC  Motor 

The  mechanical  power  developed  by  a  Brushless  DC  Motor  is  given  by  the 
equation, 

P  =  (7.062x10~3)  Twm  watts 
where  T  =  Tl+T^,  the  sum  of  the  load  torque,  T^  and  the  motor  restraining 
torque,  T^  in  oz-in  units  and  wm  is  the  motor  speed  in  rad/sec.  units.   At  no 
load  TL  =  0  and  at  stall  torque  wm  =  0,  thus  power  reaches  a  peak  value  in  the 
region  between  no  load  and  stall  torque  conditions. 

The  evaluation  of  steady-state  peak  power  can  be  obtained  analytically  in 
terms  of  the  load  torque  by  substitution  as  shown  below. 

Given  that  T=Tj__  that  is  Tl>>T^  and  T  =  ktim  with  wM=wNL  where  kt  is  the 
motor  torque  coefficient  and  kW]_  is  the  slope  of  the  speed-torque  curve. 
Then  the  power  output  is 

po  =  ktim  <w  "  kwlTL) 
and  since  steady-state  speed  conditions  are  assumed,  the  current  is  given  by 

v   -  kuw 
sg    b  m 

lm  =    RA  +  Rq 

where  vsg  is  the  voltage  at  the  terminals  of  the  power  supply,  kb  is  the  Back 
EMF  coefficient,  Ra  is  the  armature  resistance  of  2  phases  in  series  and  Rq  is 
the  sum  of  2  power  transistors  forward  'conduction  resistance.  By 
substitution, 

k  v  k  k  2 

Po  -  ^ — ^§  (w  ,  -  k  nT,  )  -  ^— -  (w  ,  -  k  t  ) 
Ra  +  Rq   nl    wl  L    RA+Rq   nl    wl  L 

and  to  find  the  conditions  for  peak  power  output,  it  is  necessary  to  compute 
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the  rate  of  change  and  set  it  equal  to  zero.   Thus, 

dP      k.k  nv       2k  kk  w  ,    2k.k.k  ,2 
o      t  wl  sg       t  b  wl  nl      t  b  wl 

dTL      RA  +  Rq       RA  +  Rq        RA  +  Rq   L  = 
Solving  for  load  torque,  TL  yields 

TL  =  wnl/kwl  "  vsg/(2kgkwl) 

TL  =  (2kbwnl  -  vsg)/(2kbkwl) 
and  substitution  into  the  power  output  equation  produces  the  conditions  for 
peak  power  output  as  follows: 

po  =  ktvsg/^Ra+Rq)[wnl-kwl(2kbwnl-vsg)/(2kbkwl)] 
"  ktkb/(Ra+Rq)[wnl-kwl(2kbwnl-vsg)/2kbkwl)]2 
After  expansion  and  cancellation  of  like  terms,  the  peak  power  output  reduces 
to 

Po  =  ktvsg/[4kb(Ra+Rq)] 

As  a  numerical  example  consider  the  following  values  of  a  typical 
Brushless  DC  Motor  operating  at  a  steady-state  speed  of  3060  RPM  (320.4 
rad/sec)   with  a  constant  terminal  voltage  of  30  volts. 

kt  =   15.9  oz-in/ampere 

kb   =   0.112   volts/rad   per   sec 

kadj   =  °-825 

Ra  =   1-37  ohms,    Rq  -   0.10  ohms 

kwl  =   1*67  rad  per  sec/oz-in 
Thus   for   this  example  motor, 

P0  -    (15. 9)(0.825)(30)2(7.062x10~3)/[U(.ll2)(0. 825)0. H7) 

PQ  ■  153.5  watts 
and  occurs  at  load  torque 

TL  =  [(320.4)/(1.67)  ~  (30. )/(2)(.112)(. 825)0.67)] 

TL  -  94.65  oz-in  =  5.92  lb-in 
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The  curve  of  power  output  versus  load  torque  for  this  example  motor  is 
shown  in  Figure  9.  It  must  be  noted  that  the  above  analysis  is  valid  for  the 
assumption  that  motor  current  produced  fields  do  not  have  a  measurable  effect 
on  the  magnitude  and  phase  of  the  flux  generated  by  the  rotating  permanent 
magnets  and  also  that  a  constant  average  flux  exists  in  the  air-gap  between 
rotor  and  stator. 
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APPENDIX    B 
CSMP  Listing 

The  system  described  by  the  schematic  diagram  in  Figure  H  is  programmed 
in  the  CSMP  language  as  shown  below.  The  computer  used  was  an  IBM  370  located 
in  the  W.  R.  Church  Computer  Center  at  the  Naval  Postgraduate  School. 
Following  the  program  listing  is  a  User's  Guide  for  the  recommended  procedure 
in  setting  the  constants  and  parameters  for  validation  with  measured  motor 
performance.  Also  included  are  typical  output  waveforms  for  a  step  input  of 
supply  voltage. 
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CN  CCE 
ICN  CC 
PICTIC 
HCTCE 
EGAD 
E  LOAD 
F  THE 
CEANIC 

NSISTC 

)     —    1 

N  S  x  S  T  U 

)  —  I 
NS I  ST 0 
E  'I  A  G  E  , 
El  AGE. 
1  B  C  J  1 1 
MO  DEE 


UE  CONSTANT 
Z-IN/AMF) 
EMF  CONSTANT 


(VCIT/RAC/S) 

MFG  SUPPEIEL 


IORQ 
1  (C 
EACK 
ANT 
FCR 
THE 
THE 
FFIC 
FFIC 
EFFI 
N  CF 
(CZ 

II 

THFC  EEEJCTICN  ^EAES 
MOTCE  SXSTEft  (SEE  EELCW, 


(CZ-IN/AftP) 

(VCET/RAD/SEC) 


{ C  H  M ) 
(HENRIES) 

ft  C  I C  hi 


MCTOF  CCNSIAN1S 
B  C  I C  F 
ft  C  T  C  F 

IFNl'cF  THE  flCICB  (CZ-IN/EAC/S) 
I  ENT  CF  THE  EC  AD  (C  Z-  1  N/  E  AD/  S) 
CIENI  CF  LCAC  THRU  EELUCT1GN  GEARS 

THE     MCTCE    SYSTEM     (SEE    B  EL  C  ri  ,  NC-S"C  R  I) 
-IN/S-S) 


it 


ii 


A  I, 

HE 

RS 

h£ 

hS 

HE 

E S  FCR  CURRENT 
VI N,  MEASURED 
VIN,  ^EASiJREC 
RESISTANCE    CE 

FC3     I  FA&SIST ICN 


NC-SCET     SECIICN) 


(SEC) 

THE  MOTCE 


TIME  CONSTANT  CF  THE  MCTCE 

ELECT  PIC AI  TIftE  CONSTANT  CF 

FCE  CUBE  EM  PATH  A-E 

EIECTpIoAL  TIME  CONSTANT  CF  THE  ftCTCB 

FCF  CUE BE  NT  PAIR  E-C 

ELECTRICAL  TIME  CONSTANT 


CF  THE  MCTCE 


E  A  I  H 

FBCfl 

FECF 


C-A 
+  TC 
-  IC 


iv1D 

tie 


ECSIIICN, 
ECSIT1CN, 


tETWEEN  5  A  I U EAT ICN  AND  CUTCFF 


B2 


ILE:  BM1BA 


csme 


A1 


IBANSI5ICR  SATURATION 
1BANSISICR  CUICFJ 


RSAT  =  EQUIVALENT  CIRCUIT  RESISTANCE  CF 
RCUT  =  EC  CIV  A LENT  CI  ECU II  fr £31 STANCE  CF 
TTIMS  —  SWITCHING  1IPE  CF  TSAKSISICB 

IHBSI  —THIS  IS  IHDEG  WHICH  IS  SET  TO  ZEEC  AT  STABI,  BESE 
AT  3o0  DEC  AND  CAN  AISC  E£  USEE  AS  SI^JIATICN  SIC 
THADV  —  THETA  ADVANCE  FOR  CCKKUTAITCN  jCC'«  IS  ECSIIIVE) 
THCCN  —  CCdMUIATION  CCNTBCL  A  NC-I  £  =  I  H  E  SI  ♦  !£>  A  E  V  IN  DEGREES 
65  INTEBNAI  RESISTANCE  CF  IKE  SUFFIX  VCIIAG5 


10  2EBC 

ANGLE. 


DSCE1 
KT=Kld*KADJ 
KB=KBM*KACJ 
BIB  =  EL/ (N**2 
JLP  =  JL/  (N**2 
J  =  JM  ♦  JLP 
E  =  fig  *•  ELF 
A1  =  LA  /  BA 
A2  =  J  /  £ 
A3  =  LA/  (BA  +RSAI) 
ES1=BS/2. 0 
BS2=ES/2.G 
VINIC  =  15.0/KINT 


2THCD  SI1FE 

(NAMIC 

^Jlf=    VINhAF 
V1B    =  -VINrJAF 
VIN  =  VIE  -  VIE 


VSGDEL 
VSGERR 
VININT 
VINHEG 
VINHAF 
VBAMP1 
VBA3E2 
VINHF  1 
VINHF2 


+3.5,  VSGlEL) 


(VSGREF  -  VSG)/2.0 
DEADS?  (  -0.5 
INTGBL  (VIMC,  VSGEaB) 
KIN  I  *  VININ1 
VINEF1  -  VINF.F2 


=  RAMP  (CO) 

=  B  A  M  P  (  j  .  0  0  1 

=  K  R  A  M  F  *  V  R  A  «  P  1 

=  KSAMF  +  VF-AME2 


) 


INTRODUCE  UNITE  TRANSITION  II 
BY  INCCFFC BATING  EXECNENIIAL  B: 
Cj#EXP  PROVIDES  A  REALISTIC  EXP 
AND  SATURATION  «£]ICH  INCLUDES 
THE  IIMI'IEB 


flE  ECB  IEANSISIOR  SWITCH 
ISE  AND  CECAY  INTO  S5i  1  —  Stt6 
ONENTIAL  IEANSIIICN  EEIBEEN  CUTOFF 
SA1UBATICN  DELAY  WhEN  EUI  T HFC  UGH 


C.1EXE 

G.2EXF 
Q3EXF 
Q4EXF 
Q5EXF 
QfaEXE 

:lerr 

5SEBR 
RQ1    = 

Bv2    = 

BQ3 

SCj4 

RCj5 

RC.6 


=     R  E  A 

=    BE  A 

=    PEA 

=    u£A 

=    RE  A 

=    BE  A 

Q1EXP= 

Q  1  £  X  P  = 

=    LIflll 

=    LIilII 

=    LIMIT 

=    LI  ill  I 

=    LIMIT 

=    LI  All 


LPL  (. 
LPL  (0 
LPL  (J 
L  F  L  ( 0 
LPL  (0 
LPL  <G 
0.  1, 

o.  1, 

{£5  AT 
(FSAT 
(  B  S  A  T 
(BSAI 
(BS  AT 
(FSAT 


.75 
.  75 
.  75 
.75 

v_Z  E 

<  E 

/  B 

$  B 

,  B 

,  B 

,  R 


TT 

,     1 

,  T 
,  T 

n 

§    a 

x?  = 

X  5  = 
CUT 
CDi 
CUT 
CUT 
CUT 
CUT 


I  d  E  - 

1 1 M  E 
TIME 
TIME 
TIME 
TIME 


LINEAR 

DIODE 

VD1D  = 

VDbD  = 

VD3D 

VD2D 

VD5D 

VD1D 

RD1 

BD2 

ED3 

ECU 

ED5 


CICDE  MODEL 
TUuN-CN  VOLTAGES 


6CC0.0 
=  faCCC. 0 

C  t  I  I  •  ■-- 

=  6CG0.0 

=  eoc 


cCOO'o  + 

EC  NSW  (VD1D 

*  C  N  S  K  (\U±Z 


VA.IND 
VCINC 

V  a  I  ft  E 

V  A I  Is  C 


VC 
VE 


F  C  N  S  4 
I C  N  S  W 

F  C  N  S  W 


VD  JE 
V  D  4  C 

VD5C 


INC 
INC  ' 

a  SAT 

BSAT 
B  S  A 1 

RSA1 
BSAT 


0.75 
0.  7 
0.7 
0.7 
0.7 
O.  / 
C3FX 
•*3EX 

♦  EC'JI 

♦  ECU  I 
+  RCUI 
+  FCUI 

♦  ECU! 
*PCU1 


-SW 
5-S 

C  _  Q 

5-S 

5-S 
5-S 

F=G 

P  =  C 

*C1 
*C2 

*Q2 

*£4 
*v5 
*Cc 


A  J 

W  *4  * 

A  ' 

.  1 

.1, 

EXt 
EXP 

EX 
EXP' 
EXP' 
EX? 


VCINC 
VEJNC 

V  A  I  N  C 
VCINC 

V  E  I  N  0 

V  A  I  N  C 
ECUT# 
fi  C  'J  1 1 

C  /"«  P  fT 

i*.  C  L  j.  f 

R  CUT, 

RCUT, 


BCUT 
FCt  I 
ECU  I 
BCCT 

FCLi 


kl4EXI  =  0.1,     C.5EXF=C1,    Cfc£XE  =  0.1 
C4EXP  =  C.1/     Q5EXF  =  C.I,     ^cEXE  =  C1 
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FILE:     BHIBA 


L  S  -'•!  F 


A1 


BDo     =    fCNSfc(VLi.Cf     ESA1,     BCU1,     FCII) 


*  THE  1 

*  VGLIA 

*  APPLI 

*  THEVE 
V  A  N  T  B 
VBNTH 
VCNTfl 


HEV3NIN  EQUIVALENT  CISCUIT  METHOD  IS 
G£  ['ID-FCINl,  N,  A3  EErEEEKCE  ANl  THE 
EC  IC  COMPUTE  THE  PHASE  CIFFEMS. 


USED 

LCCr 


Vi  I T  ri  I  fi  E  S 

MEIriCC  CF 


UFFLY 
ANALYSIS 


MIN  ECUIVAIEM 

=  VNl*FEvAS2/ 
=  VN1*E£C.£S2/ 
=  VN1*FEC.C32/ 


VCITAGES 

FE£A5l  +  FEC.ASfc)  +  V  N  2  *  a  E  C  AS  1/  (  BE  C  AS  1  ♦ 

EE£BS1*FE£ES2)  +  V  N2*  F  E  £  ES  1/  (EECES  1  ♦ 

6EQCS1+REC.CS2)  *  VK2*EE£CS  1/  JBEC.CS  14 


VN2*BECAS  1/  (EECAS1  +  B 

\,  k    )  *  t  t  r  E  ci/Jcr/-t.ci  +  c 


*    TEANSISICr/DJLCE    E'.UIVriLENT     RESISTANCE 


REQA1 
REQA2 
RE0.E1 
EEQB2 
REQC1 
BE£C2 


F£1 
F^ 
F£3 
FQ4 
h£S 
FK6 


FC  1/  (EQ  1+RC1) 

B  E  2  /  (  F  £  2  ♦  F  D  2  ) 

FC J/  (F£2+ -D j" 
ECU/  (KM  +  F.D4 
BE 5/  (F£5+RD5 
r*E6/  (Fw6  +  BD6' 


COMPUTATION  OF  IHEVENIN  E£UIVAIENI  FESISTANCE 


EEQAS1  = 

BEJ.AS2  = 

BEO.BS1  = 

EEQBS2  = 
BE£CS  1 
EEQCS2 


FE£A  1 
B  E  Q  A  2 
B  E  «  E  1 
SE£c2 
BEQC1 
FEQC2 


IHEVENIN  E£U 
BEQA  =  FE2A" 
BE0.B  =  BECE- 
BEQC  =  FE£C31 


FS  1 
rS2 
FS1 
FS2 
BS1 
BS2 


gc^  +  RSz 

UIVALEK1  RESISTANCE 
S1  *  FE£A^2/(RE£AS1+FL£AS2) 
S1  *  FEvfc£2/  (FEwBS  1+FE£ES2 
Si  *  EEQCS2/ (EECCS1+EECCS2) 


TOTAL  LCCE  BEblSTANCE 

REQAE  =  2.+BA  +  EEC.  A  ♦  FEQE 

BEO.BC  =  2.*BA  ♦  BECE  +  BEQC 

BE£CA  =  2.*Rii  ♦  EEQC  ♦  RECJ 

LOOP  CUSFENT-  MAXIMUM  VALUE 

IAS  =  (VANTH  -  VENIE)/£E0.A£ 

I  ABA  =  »(VI^  +  VEMED  -  VEM1-A)/(FS  ♦  2.*  3  A  ♦  BEGjB2 

IBC  =  (VcMK  -  »CNTE)/FE£BC 

IBCA  =  -{VIN  +  VEKFB  -  VEMEC)/JFS  +  2.*FA  +  BEQC1 

ICA  =  (VCN1H  -  V  ANTE) /EEQC A 

ICAA  =  -(VIN  +  VE1FC  -  VEMEA)/(SS  +  2.*FA  ♦  FEQA1 


+  FE£  A  1) 
+  FE£E2) 
♦  FE£C2) 


TIME  CONSTANTS 

AEIAU  =  ^.*  LA/BE£AE 

BCTAU  =  2.*    LA/EEC  EC 


CATAU  = 


*  LA/EE£CA 


0,AETAU,1AL 
0,ECIAU,IcC' 
0,CAIAU,ICA 
"  E  E  M  S 


LOOP  CUFFENTS 

IMAB  =  SEALPL  (0 

IMBC  =  BEftLPI  C 

I MCA  =  BEALPL  {0 

NET  LEO  (PHASE)  C  0  u 

IMA  =  IK  At  -  IMCA 

1MB  =  Ifldv.  -  IMAb 

IMC  =  1KCA  -  IMEC 

IMCD  =  DEFIV (0.0, IMC) 

IMCDB  =  TCIMC+IMCE 

IMCAC  =  hEALPL  (0. 0, TCIMC 

IMF  =  FEAL?L(0.0,G 

VSGF  =  EEALFL  (J. D, 


IMCDE) 
0.C05, VSC) 


INDUC1CB  ljSUCEL  VCL'TA: 
CIADT  =  DEiilV  <C.C#  l?.b) 
DIBDT  =  DEFIV  (CO,  I«E) 
DICDI  =  DEB  IV  (C. 0,  IMC) 
VAIN£  =  LA*I  IAD  I 
VB1NC=L A*C IEDT 
VCINC  =  LA*CICD  I 
PHASE  CHM«S  LAy  VCLIAGE 


ES 


B4 


L  E :  E  a  1  &  a 


CSMP 


Al 


VAOLV=IMA*  <&A 
VEOLV=IKE*  (BA 
VCOLV=IKC*  (BA 
PHASE  VCLlaGE 
VAO  =  ISA*  (BA 
VBO  =  IMB*{RA 
VCO  =  IMC*(RA 


♦  LA*C1ACH-VEMFA 

♦  LA  +  CIECH-VEf.FE 
♦LA*EICEI+VE£FC 


CCMPUT 
WHERE 
VSG  = 
VAG  = 
VEG  = 
VCG  = 
VAC  = 
VOLTAG 
VQD1  = 
V^D2  = 
VG.D3  = 
VCD4  = 
VQD5  = 
VQD6  = 

PVAG  = 
PVBG  = 
PVCG  = 
PVAC  = 
PVBA  = 
PVCB  = 


AUG 

VIN 

VIN 

VSG 

VSG 

VSG 

VAG 

IS  A 
VSG 
VAG 
VSG 
VEG 
VSG 
VCG 


N  OF  PCWEB  SUPPLY 
IS  SUPPLY  VCLTAGE 

-  IK+ES 

*  P.  EC  A  2/  {FE£A1  ♦ 

*  R5^E2/  (5ECE1  ♦ 

*  EECC2/(E£C.C1  ♦ 

-  VCG 

CROSS  1EAN5ISICSS 

-  VAG 

-  VEG 

-  VCG 


1  E  5  ft  I N  A  I 
ANC  1ft     IS 


B£QA2 
FECE2 
a  E  £  C  2 


ANE  EICCES 


VCLIAGrS  (GND  BEE) 
CUH" 


b  U  F  1 1  Y 


BESI 


LIMIT  1-60. , 

LIMII   -63-, 
LIMII  1-bQ., 
EVAG-PVCG 
EVEG-PVAG 
FVCG-PVS3 


ec 
ec 


id 

IQ2 

IQ3 
IQ4 

IQ5  = 

IC6  = 

pqi  = 

FQ2  = 

PQ3  = 

pga  = 

FQ5  = 

PO  = 
PQTOT 


FCNSr* 
r  C  N  S  * 
iCNSb  i. 
ECNSw  (£ 
CNSW  (< 


LV.IIJB   1  J  «  ^  , 

=  FCNSki  (S»3, 

cr":- '  s  w  a , 

SW5, 
{ S  *'  6  . 
* 


F  C  N  S  »' 

VEC  1 
VL-C.2 

VEQ3 

VDC^ 
VEQ5 

VC^b 
=  PQ1 


IQ1 
IC3 

IC5 

Ivb 

♦  EC2 


.003, 

.003, 
.CC3, 
.003, 
.003, 

1  ""  "5 


VAG 
VEG 
VCG 


.002, 
.002, 
.303, 
.003, 
.003, 
.003, 


PQ 3  ♦  EC4  +  EC 


FQe 


iEMFb  =( 
2EMFC  ={ 
YCRMALI2 


THE  PJAS 
BEMFA  =  ( 
EEMFb  = 
E 
N 

VEMFA  = 
VEMFE  = 
VEMFC  = 
VEMFAC  = 
VEMFE  A  = 
VEMFCE  = 


A  N  G  L 
*S1N 
♦  SIN 
*SIN 
LEG 
BEMFA 
fcEMFE 
EEMFC 
VEME  A' 
VEMFE- 
VEMFC- 


LCCK 
.♦TEE 

.  *1HE 

.*lhE 

E 

KB/.KK 

KE/KK 

KB/KK 

EM  EC 

EMEA 

E*EE 


CCNFUSEC  BU: 
LA*  (11*Pl/t) ) 
IA+  (7*?I/b" 
I  A*  3*PI/6 


IHE  PHASE  BEIAIICNSEIf  IS 
♦.59*SIN{1C.*lHElA*(1*Fl/6 
))*.5^*SIJi<1" 
J5+.5S 


A*  I  1 

♦  (5* 

9*SIN  (10.*TdEIA*  (9* PI/ 


5y*SIN  (10.*THE1A*  (5*£1/ 


!ii 


PIG  HI 
) 


3  *  inK 

3   *   rHLV 

3  *  fof 


1CRQUE  IC  MAKE  MC1CF  1URN 

TA  =  IMA  *  Kl  *  EEMFA/KK3 

IB  =  IME  *  KT  ♦  EEMFE/KK3 

TC  =  IMC  *  Kl  *  EEMFC/KK3 

TM  =  1A  ♦  Tfi  +  IC 

1MM  =  Kl  *  I M 

IBM  =  EM  *  WM 

VKBfcfi  =  KE*WK 

IORCCE  IN  LEGE  A L  VALUE 

1MINI=INIGBL (0. :,1M) 

1L  =  1LL+SIEP  (1.0E-3) 

IN  2  =  TNI  *  M.O/t) 
HM  =  &EALEL  (C.0,A2 


UMINT 


,  1  ht  2 
I  NT G i? I  (  J  .  0  ,  »  ft  ) 
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file:  E  :i  I  r-  A 


LS.U 


A1 


WB=32G 

WfifcPB 

HflRPBE 

THEIA 

IHDEG  =  IHETA 


.  a  a  z  5 

=  iiK  *  1 3  0  .  /  P I ) 

=  KBRPB/N 
=  IMGfcL(O.O.hfl) 
*  (18C.G 


/PI) 


COMMCTriilCN  ADVANCE  HEBE 
THCON  =  1HRSI+IHADV 

PWRMI  =  iJB  *  IB  *  .3C7D625 
P  W  R  M  =  E  B  *  .0070675  *  tf  B  *  *  2 
P  W  R  =  EI  *  . 0 0 7 0 c  7  5  *  »  B  *  * 2 
II  =  VIN/  (2*8X01) 

THIS  PROCEDURE  PROVIDES  A  SIBPIE 
DIRECTION  -JHIIE  KEEPING  THE  LCAC 


BECHANISB  EOF  EEVEESING 
AS  AN  CEECSING  1CEC.UE. 


1HE  BC1CB«S 


PROCEDURE  INl  =  FWCEfaC  ( fc  B , T  B , II) 
If jlfl.IT.G.O)  GO  IC  1 

T  N  1  =  IB  -  II 

GC  TC  b 
1   III  =  IB  ♦  II 
5   CONTINUE 
ENDPRCCEDUBF 

*  IHIS  P50CEDUBE  EESE3S  THE  V  A E I J E I E  1HRSI  IC  C  AFTER  EVEEY  360  CEGEEFS 

*  OF  MECHANICAL  ROTATION-   1EIS  Ii  FUNDABEbilAI  TC  THE  SIBLLA1I0N  OF  All 

*  SWITCHING  AND  EtSIlICN  SENSING  AC1ICN. 


PROCEDURE  1HRS1=B 
1HES1  =  ABCC  ( 
IF  (I  HE  ST  .11. 

ENDERCCECUEE 


ESET(IEDEG) 

IFDEo,  joJ.) 
0.0)  TiiRSI  =  In  EST  *    360. 


*  1H 

*  SE 

*  RE 
PROC 

IF 
GC 

4b 

46  C 
IF 
IF 
IF 
IF 
IF 
IF 
10 


11 

12 
13 

14 

1b 

20 
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NSCR 

GARD 

EDUE 

(THC 
TO 

IHCO 

ON1I 
THu 
THC 
THC 
I  EC 
THC 
THC 

SE1 

S  £2 

SE3 

GOIO 

SE1 

SE2 

SE3 

GCIO 

SE  1 

SE2 

SE3 

GCIC 

SE1 

SE2 

SE3 

GOIC 

SE1 

SE2 

SE3 

GCIO 

SE1 

SE* 

SE3 

CONI 


FCC 

S  S 

L  ES 

E  S 

ON 

46 

N  = 

NUE 

ON 

ON 

CN 

ON 

CN 


E'JUKE  WAS  IAKEN  FEOfl  Ihl  NEXI 
HOUID  GE NEBAIE  A  SINGLE  VAIOF. 
S  OF  DIRECTION  CF  RC'IiUCN 
fc1. SE 2. SE3  =  H AIL  (IHCCN) 
.GE.130.)  GC  10  4  5 


EOCEDURE 
f  UNO  HON 


EECAUSE  IKE  HAII 
CF  ECS1IICN 


IHCCN  -  18u 


GE. 

GE  . 
GE. 
GE. 

G5. 


33. 
60. 
9  0. 
20. 


AND. 
AND. 
AND. 
AND. 
AND. 
AND. 


IHCCN 
THCCN 
IHCCN 

IHCCN 
IHCCN 
i  H  C  0  N 


IT.  J 


IT. 
II. 

II. 
II. 
II. 


9  0. 
1*0. 
150. 

1d0. 


GC 
GC 
bC 
GC 
GC 
GC 


IC 
10 
IC 
IC 

11 

IC 


10 

1 1 

12 
13 
14 
15 
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2C 
=  1. 
=  0. 
=  0 . 

2C 
=  1. 
=  1. 
=  0. 

*0 

=  c. 
=   1. 

=  0. 
20 

-  VJ . 

=   1. 
=  1. 

2  0 

-  0 . 

=  c. 
=   1. 
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IE:  EMIRA 


csa, 


A1 


DERCCECURE 


THIS  P 

DETEEM 

THE  VA 

APPRCX 

VCITAG 

SMI  1H 

SIMULA 

THCON 

IS  I«P 

BEGINS 


F.QCEDURE  SIMULATES  GfrNERAI  C  Cf.M  UT  A  1 IC  N  , 

IN  ES  IHE  AIGEEBAIC  SUM  CF 

FIAELE  ELUX  CCMfCNENIS  FC5  USE    IN  CCHEU1IMG  THE  MCICF'S 

IflATE  BACK  EtfF  (BEMFT).  II  AISC  SUETBACIS  THE  GEMEAIED 

E  FBCM  THE  PBCEEB  SUIPI*  VCITAGF. 

EU  SW6 

TE  POWER  TRANSISTOR  1BIGGE3S  EEING  ENEEGIZED  CR  SfcllCHED  OFF 

IS  THE  VAEIAELE  THFCUGH  hHICH  1HE  SWIICEIKG  ICGIC 

LEMENTED.  REVTIM  IS  IKE  IlrfE  AT  WKICE  CICCKWISE  COMMUTATION 


OCEDUB 
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IFJTIM 
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IF  (T 
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IF  THC 
IF  (THC 
IF  (THC 


C 


81 

SW2 

SW5 

SW4 

SW5 

SW6 

BEME 

VN1 

VN2 

IM  =  I 

GC  I 

SU1 

SW2 

SW3 

SWM 

SH5 

SW6 

BEMF 

VN1 

VN2 

ia=i 

GC  T 
SW  1 
SW2 

SH3 

sua 

3W5 

S'it 

BEME 

VN1 

VN2 

IM  =  I 

GO  I 

SK1 

SW2 

SW3 

swa 

SW5 

swe 

BEfiF 


E  SS  1 

HCCN, 
E  .GT 

in  = 

HCON  1 

CN1  . 

CN1  . 

C  N  1  . 

GN1  . 

CiN  1  . 

CN1  . 

ISE 

INUE 

CN1  . 

CN1  . 

CN1  . 

ON  1  . 

GN1  . 

CN1  . 

=    C. 

=  0. 

=  G. 

=  1. 

*  1. 
0. 
=  B 
VIF 
VIE 


,3A2 
BEKF 
.REV 
AMOC 
.IT 
GE. 
GE. 
GE. 
GE. 
GE 
GI 


1 

COM  (J 

.GE. 
.GE. 
•  Gil. 
.  GE. 
.GE.  1 
.  GE.  1 


,SfcJ, 
A,EEM 
II  fl  ) 

{THC 
.  C. 
0.. 
30.. 
bG.. 
90.. 
20.. 


SW 

FE 

G 


C 
A 
A 
A 
A 
A 
u.  .  A 
T  A  1 1 C 


ND 
NE 
ND 
a  I 
ND 
NC 
N 


4,S 

,BE 
C  1 
18 
THC 
.  T 


W 


IH6.EE  KFT 


0. 
30. 

eo. 

90. 

2C. 


A*C 
AND 
AND 

AND 

AND 
AND 


I  = 


E^EC 

-  V 

-  V 


-  E  E  M  F  t 
EMFC 
EMEE 


ttC 

c 


MA 
C 


fcC 

=  1. 

=  c. 

=  0. 

=   1. 

=  c. 

=  0. 

1  =  EEMFfl  -  EEMFE 
=  VIF  -  VEMEA 
VIE  -  VEMFE 

6C 
1. 

0. 

c. 

0. 

G. 

1. 

=  EEMFA  -  EEKEC 

VIF  -  VEMFA 

VIB  -  VEMFC 

cG 

=  G. 

=  0. 

=  1. 

=  C. 

=  C. 

=  1. 

T  =  EEMFE  -  EEtfFC 


MA 
C 


MFC. 

0  10 
0.) 

0N1  =  THCCN1  4  180. 

HCON1  .IT.  3G.)  GC 

HCCN1  .IT.  60.)  GC 

HCON1  .IT.  90.  GC 

HCCN1  .11.120.)  GC 

HCCN1  .11.150.)  GC 

HCCN1  .11.18C.)  GC 


VN  1  ,VN2    =  CCMM  illlE, 


VEM£A,VEMFE,VEMFC) 
0 


IC 

IC 
IC 
IC 
TC 
TC 


THCON  1 
IHCON  1 
THCON  1 
THCON1 
THCON  1 
THCC-N1 


IT 
L'i 

IT 
II 
II 
II 


30. 

fcO. 

90. 
120. 
150. 
130. 


GC 
GC 
GC 
GC 
GC 
GC 


TC 
TC 
TC 
TC 

TC 
IC 


50 
51 


5a 

55 


5a 

C  C 

50 
51 
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FILE:  E.1T&A 


CSMF 


A  I 


5a 


VN1  =  VIF 
VN2  =  Vic 
IH-INE 

GC  1C  bO 

SW  1 


-  V £  M  E  E 

-  VErfEC 


55 


SW2  = 
SS  3  = 

swa  = 

SW'5  = 
5W6  = 

B  E  «  F  I 
VN  1  = 
VN2 
IM=IME 
GO  10 
SW  1  = 


u  • 
1. 
1. 

0. 

c. 

0. 

=  BErtFE  -  EEMEA 
VIF  -  VEMrS 
-  VEMFA 


-     VIE 


cO 
0. 
1. 
G. 
0. 

1. 

0. 


60 


SWz  = 
Sii3    = 

swa   = 

SW5  - 
sut   = 

BEflFT 
VN1  = 
VN2  = 
IM=IMC 

CONIINUE 


=  EEHFC  -  EEMEA 
VIF  -  VEMFC 
VIE  -  VEMFA 


ENDPBOCEEUEE 


TERMINAL 
TITLE 

TIMEE 

FINISH 
PRINT 
VAIND 
VSG,  V 
IAB,I 
1.1  A. 
LABEL 
OOIPUT 
PAGE  X 

E 
BESET 
LABEL 
CUTFUT 
PAGE  X 

£ 
LABEI 
OUTPUT 
PAGE  X 

E 
LABEL 
C  U  T  P  U 1 
PAGE  X 

E 
LABEL 
OUTPUT 
PAGE  X 

E 
LABEL 
CUIPUT 
PAGE  X 

E 
LABEL 
OUTPUT 
PAGE  X 

E 
LABEL 
OUTPUT 
PAGE  X 

fc 
LAEEL 


BASIC 

FINIIM 

CEIMIN 

T  H  E  S I 
*  A  .  W  tl 
,  V  fc  li* 
AG, VEG 
AEA, IB 
1Kb,  I 

EBASE 

TIME, 
YELCT 
NC 
ERINT 

PC  WEE 

TIME, 
YELCT 
NC 
/ICIb  S 

1 1  ri  E  , 
YFLCI 
NC 

POWER 

THE, 
v  P  L  C  i 
NC 
TOTAL 

TIME, 
YE  LCI 
NC 
EE-rIL 

1 1  fl  E  # 
YELCI 
NC 
VOLTAG 

TIME, 
YFLC  I 
NC 

LP-FI 

TIME, 
Y  E  L  C  I 
NC 
cACK  E 


DC  MOTCE  SYSTEM 

=    .050,  CITDEL  =  .QC00275,  ERCEL  =  .00C275, 

=1.3E-10 

=    40J.C 
hFM,     IHESI/IM.IMINT.IMINT,  *  MI  NT  £  I  A  V  ,  I  A  V  ,  W  MA  V  ,     .. 
D,     VCIND,VN  1. VN2,     V  AC , V  EC , VCC  ,  C J , V QD H , V Q D 5 , V £C 6 , 
,  VCG,VDlD,VB2D,VD5C,VC4i;,VC5£,VC6C,     ... 
C,  IECA,1CA,ICAA,IM,IECEG,     T HCC N  ,  IhCO N  1 ,     ... 
MC,     E*E,VIN,VINHAF,IMCAC,VSGF,TMM,TEM,VKEWM 

C  CUE  BENT  -  AMES 

IMC,SW2 


SUPPLY  CURRENT  TC  MCIOE  -  AMPS 
I  M  ,  S  *  4 


SEEZC  -  BE* 
W  P.  R  t  M  ,  £  *  5 


S  U  P  P I i    TERMINAL  VCITAGE  -  VCLiS 
VSG,£  /it 


DEVELOPED  IUEQUE  -  CZ-IN 
T  M  ,  S  W  1 


TEBEC  PEASE  C  lURREM  -  A:1FG 
IMCAC,SUi 


E  ACROSS  A-C  TERMINALS  -  VCLTS 
VAC,SW3 


L TEBEC  PCwOR  SUPPLY  CURRENT  -AMPS 
I  fi  f  ,  S  fc  5 


:1E  F  Li  AGE  C  -  VCLTS 
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LE:  BMTEA 


CSM? 


A1 


-  HAT1S 


3UTPUI  1IHE,  V£HFC,SW  1 
PAGE  XiPICT 

END 
LABEL  Ifl  CLVELCEED  TCBsJUE  - 
3UTPUT  LIME,  1{1fl,SW3 
PAGE  XY£LC1 

END 
-ABEL    ECBEB  CUTCU1 
IfllPUT  1IKE,  P«R#SW5 
PAGE  XYflCl 

END 
.ABEL  LP-FILTEFED 
)UTPUT  1IHE,  VSGF 
PAGE  XYtLCl 

END 

SICE 
DJOB 


C2-IN 


SUPPLY 
,SM1 


TEfiMINAL  VCL1AGE  -  VCL1S 


)ISSFCP.MYDInS 
DRAii=  1-ENE, 
MODIFY=1-£NE (J 
«QDIFY=2-ENC (3 
MODIEY=3-END  <J 


C  * 


(SIZE  =  8.  -  j.)  , 

(CV££FLCi,CCt 


KNEE=0. ,3..SIZE=8 


,3. 
(OV£FELOl,CCfiliEE=0.,3.,SI2£  =  8.,3 


:}• 
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Users  Guide 

The  following  procedures  are  recommended  to  be  used  in  running  the  CSMP 
language  program  given  in  Appendix  B: 

1.  Obtain  the  motor,  load  and  power  supply/conditioner  data  from  the 
manufacturer  and  enter  the  values  in  INITIAL  section  at  the  beginning 
of  the  program. 

2.  Set  parameter  KRAMP  to  the  desired  input  voltage  V  as  follows: 

V    3 
example  V  =  30  volts,  KRAMP  =  -  x10  =  15000.   Under  these  conditions, 

a  fast  terminated-ramp  input  results.  This  is  equivalent  to  a  step  input 
of  V  volts. 

3.  The  constant  KADJ  must  be  set  experimentally  to  produce  the  measured  motor 
Back  EMF  voltage.  To  do  this,  remove  the  asterisk  in  the  statement 

*  WM  =  320.4475 
where  320.4475  is  the  magnitude  of  the  driver  motor  speed  in  rad/sec.  (in 
this  case,  speed  =  3060  RPM).   Run  the  program  with  different  values  of 
KADJ  until  the  Back  EMF  across  2  windings  agree  in  magnitude  with  the 
measured  value  obtained  from  the  motor. 

4.  Set  N  to  the  value  of  the  speed  reduction  from  motor  shaft  to  load 
shaft. 

5.  Set  RS  to  the  value  of  the  power  supply  interval  resistance  and  LA  and 
RA  to  the  per  phase  winding  values. 

6.  Set  RSAT  and  VSAT  values  for  the  power  transistor  "ON"  values  and  RCUT  to 
the  cut-off  resistance  value. 


BIO 


7.  The  value  of  BM  should  be  calculated  from  the  no-load  data  as  a 

preliminary  value  and  adjusted  experimentally  until  no-load  operation  of 
the  model  agrees  with  measured  no-load  values  for  the  motor.  Compute  BM 
as  follows: 

BM  =  (KT)(KADJ)(IMNL)/(WMNL) 
For  example 

BM  =  (15.89)(.63)(0.3)/(320.i4) 

BM  =  0.00937  oz-in/rad  per  sec. 

Notes:  1.  The  asterisks  in  column  1  of  the  program  statement  makes  the 
statement  a  comment  and  disregarded  by  the  CSMP  Translator.  These  statements 
are  the  result  of  program  development  and  are  included  for  future  application 
in  the  development  of  the  complete  actuator  model. 

2.  A  set  of  typical  output  waveforms  are  included  here  to  indicate 
the  plot  output  capabilities  of  the  program  which  are  in  the  TERMINAL  section 
of  the  CSMP  listing. 

3.  The  load  applied  to  the  motor  should  be  adjusted  over  the  range  of 
no-load  to  near  peak  load  and  the  performance  curves  of  Figure  12  plotted. 
The  CSMP  simulation  output  can  then  be  checked  against  these  results  by 
entering  BL  as  follows: 

BL  =  TL/WM  oz-in/rad  per  sec 
where  TL   is  the  load  torque   (oz-in)   WM  is  the  motor  speed    (rad/sec). 


Bll 


POWER  OUTOUT  -  WATTS 


i 


n  =  PWR 
o  =  S!¥5 


LEGEND 


0.000     O.OOB         0.010         0.018 


0.020         0.028         0.030         0.038 

TIME 


0.040         0.040         0  080 


TOTAL  DEVELOPED  TORQUE  -  OZ-IN 


(V/WYft  N 


Ps/V 


a  =  TM 

o  =  SWl 


0.000     O.OOB         0.010         0.01S         0.020         0.02S         0.030 

TTME 
MOTOR  SPEED  -  RPM 


0.030         0.040         0.040         0.080 


LEGEND 


a=WMRPM 


LEGEND 


O.OOO     0.0O6         0.010         0.015         0.020         0.026         0.030         0.036         0.040         0.048         0.06O 

TIME 

Figure  Bl  Step  Input  Response  Waveforms 
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LP-FILTERED  POWER  SUPPLY  CURRENT  -AMPS 


a=  IMF 

o  =  SW5 


LEGEND 


0.000     0.006         0.010         0.010 


0.020         0.028         0.030         0.030 

TTME 


0.040         0.040  0.000 


POWER  SUPPLY  CURRENT  TO  MOTOR  -  AMPS 


a  =  IM 


0.000     0.000         0.010         0.010  0.020         0.020         0.030         0.035         0.040         0.040         0.000 

TIME 

PHASE  C  CURRENT  -  AMPS 


LEGEND 


s 


MM 


Ha 


MM 


HaJ 


os  DkC 


0.000     0.000         0.010         0.010         0.020         0.020         0.030         0.035         0.040         0.040  0.060 

TIME 

Figure  B2  Additional  Step  Response  Waveforms 
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